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Top partner searches are one of the key aspects of new physics analyses at the LHC. We correct 
an earlier statement that supersymmetric top searches based on decays to semileptonic tops are not 
promising. Reconstructing the direction of the boosted leptonic top quark and correlating it with 
the measured missing transverse energy vector allows us to reduce the top pair background to an 
easily manageable level. In addition, reconstructing the full momentum of the leptonic top quark 
determines the stop mass based on an Mti endpoint. 

Searches for top squarks at hadron colliders aim at a fundamental questions of electroweak symmetry break- 
ing: if the Higgs boson should be a fundamental scalar, how can its mass be stabilized? In particular, is the 
Higgs mass protected by some symmetry? Such symmetries typically predict the existence of a top partner, like 
in supersymmetric or little Higgs models [TJ [5] ■ In such a case, studying the properties of top partners allows 
us to unravel the nature of such an underlying fundamental symmetry protecting the fundamental Higgs mass 
at the scale of electroweak symmetry breaking. 

At the Tevatron, low-mass stop searches look for loop-induced stop decays [3] to charm quarks and the 
lightest neutralino 0]. Slightly larger stop masses more in agreement with theoretical expectations makes it 
more promising to look for decays to a bottom jet and the lightest chargino [5]. For leptonic chargino decays this 
signatures is irreducible from a leptonic top decay. Finally, if the stop becomes heavier and the strong decay in 
a top quark and a gluino is not yet kinematically allowed, the stop can decay into a top quark and the lightest 
neutralino [3J. Because it includes on-shell top quarks as stop decay products this channel is the hardest to 
extract at the LHC. For hadronically decaying top quarks we know that this final state has the advantage that 
we can fully reconstruct its kinematics, which puts us into a promising position to study angular correlations 
in the stop pair final state |B]. Extending this analysis strategy to semileptonic stop decays is the aim of this 
study. 

There have been several suggestions as to what we might be able to say about the nature of the stop based 
on a momentum reconstruction of its visible decay products [7]; the first ingredient to any such analysis is a 
way to extract stop pairs from the overwhelming top pair background, the second the reconstruction of the 
stop 4-momentum. Over almost the entire supersymmetric parameter space standard analysis methods fail at 
both tasks. On the other hand, fully hadronic top pairs we can tag and fully reconstruct using fat jet analysis 
tools [6, 8, 9J. Using the HEPTopTagger [6j |9] we can extract the stop signature for a wide range of stop 
masses, reconstruct the top 4-momentum with high precision, and reconstruct the stop mass using a simple 
Mt2 endpoint analysis [5] . Obviously, we can use the same strategy for the hadronic side of signatures related 
to semileptonic top pairs. 

The key idea behind the reconstruction of hadronically decaying heavy particles using fat jets is to study the 
clustering history of a jet algorithm for geometrically large jets including boosted heavy particles. This idea 
was originally developed for W bosons and top quarks [101 E] and very successfully applied to Higgs decays 
to bottoms [12], hadronic top decays [6j 8, 9! as well as supersymmetric particle decays [13]. 

The obvious question is what we can learn from boosted hadronically decaying top analyses for the equally 
boosted leptonic top decay [H] . The main shortcoming of the traditional W and top mass constrained recon- 
struction is that we cannot use the missing transverse momentum vector as a cut against backgrounds after 
we use it for the reconstruction. The aim is to instead limit ourselves to the particular features of collinear 
leptonic top decays to reconstruct the top kinematics [TJl [T5] . 

In this paper we will proceed in two steps: first, we show how to reconstruct the direction of the top 3- 
momentum exploiting the typical kinematic features of boosted leptonic top decays. A cut on this top direction 
serves as the main ingredient to the extraction of semileptonic stop pairs. Second, again using typical features of 
boosted top decays we reconstruct the full 4-momentum of the top quark. In analogy to hadronically decaying 
tops this reconstruction could be a starting point to many more phenomenological top pair analyses. 
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I. SEMILEPTONIC STOP DECAYS 

To illustrate the power of reconstructing leptonic top pairs we show how we can extract the stop signature 

pp -> Ul -> (tx?) -> te?) + (w-^x?) (i) 

including missing energy from top decay and from the dark matter agent from top pair production. Our analysis 
will proceed based on two strategies: first, we suggest a simple approximation of the top decay kinematics to 
reconstruct the direction of the top 3- momentum (Section |n| and the full 4-momentum (Section |III) . This 
way we arrive at a two-step strategy to distinguish the stop signal from top pair backgrounds based on the 
reconstructed top direction and on its reconstructed 4-momentum. 

Typical semileptonic analyses start from a set of acceptance cuts requiring at least four jets and a charged 
lepton [6 . In our case we include a hadronic top tag, which means we require 

1. exactly one trigger lepton (px > 20 GeV, \r/\ < 2.5 and isolation ExjEx t < 0.1 within R = 0.2) 

2. missing transverse momentum j>x > 150 GeV 

3. one tagged hadronic top (mass drop, p T > 200 GeV, HEPTopTagger default [5J, Fastjet [17] ) 

4. one b tag among the leading 3 jets outside the tagged top (C/A with R — 0.5, pt > 25 GeV, \r]\ < 2.5) 

5. a bottom-lepton invariant mass allowed by the top decay kinematics < \J — = 154.6 GeV. 

The missing energy cut in combination with the lepton requirement safely removes all QCD backgrounds, 
including b pair production. In our analysis we do not take into account fakes from pure QCD events. 

The main Standard Model background after these first analysis steps is IF+jets production, with an inclusive 
W cross section roughly 200 times the top pair production rate. The probability to fake a tagged top from 
W+jets events ranges around 4%, after requiring to observe a fat jet with px > 200 GeV [6j. The tagging 
efficiency and mis-tagging rate for the HEPTopTagger have to be confirmed by a full experimental study 
and should in addition be tested on a sufficiently large tt sample at the LHC. Including the additional b tag 
we expect semileptonic ti pairs to become the main background to our stop searches. The bottom jet isolation 
from the top fat jet is defined geometrically. For the b tag we assume a (mis-) tagging probability of (60%, 2%) 
as an event weight factor. Charm jets with a generically higher mis-tagging probability should not spoil our 
analysis because of their suppression compared to high-multiplicity gluon production. 

The status after these first five analysis steps we show in Table [I] The stop with different masses we assume 
to decay with a 100% branching ratio to a top quark and a 98 GeV lightest neutralino. The NLO normalization 
of the Herwig+4- [IB] signal sample we obtain from Prospino [15]. The top pair sample from Alpgen- 
Pythia we normalize to the approximate NNLO rate around 918 pb [2U]- For subleading IF+jets background 
we use the default Alpgen-Pythia normalization [21]. For signal and background we include initial and final 
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Table I: Signal (for different stop masses) and backgrounds for the base analysis. All numbers given in fb. For VK+jets, 
we only show numbers after pre-selection (rijets > 3). 
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state radiation, hadronization and underlying event. The top and stop samples we generate inclusively without 
restricted decays. 

Table [T] shows how the main analysis steps act on the signal and the dominant backgrounds. Indeed, we 
see that while the lepton and missing energy requirements are easily passed by all backgrounds, the hadronic 
top tag and the b tag significantly reduce the W+jets backgrounds and leave us with the dominant top pair 
contamination. For the hadronic top tag we observe the expected large efficiency for heavier stop samples as 
compared to continuum tt production, which comes from the harder px,t spectrum. 

While these initial results are promising, Table [I] clearly indicates that we definitely need another set of 
analysis steps to improve the signal-to-background ratio to a manageable level, to be able to cope with systematic 
experimental and theoretical uncertainties. 



II. RECONSTRUCTING THE TOP DIRECTION 



The main difference between the reconstruction of a hadronically decaying top and a leptonic top decays is 
the significantly smaller number of observables. Aside from the boost of the top and its correlations with other 
particles in the final state we have to rely on three measurements 

Ei E b mu , (2) 

assuming that we know the bottom and lepton 4- vectors after identifying the particles. The invariant mass can 
be replaced by the angle between the lepton and b momenta. 

There have been a few attempts to find strategies to identify boosted tops decaying leptonically. One of them 
are 'stuck leptons', referring to leptons which are not well isolated from the & jet from the top decay [TSIIH]. The 
invariant mass of the visible top decay products mu is not necessarily close to its upper bound if the neutrino 
is not observed. Nevertheless, the mass drop m\jm\i ranges around zero for boosted tops and VF+jets while 
it tends towards one for heavy flavor jets. Because Eq.|2| tells us that there exist three distinctive kinematic 
observables we can combine the mass drop measurement with preferably sizeable Ai?^ and Ee/(Eg + E b ) ~ 0.4 
to construct a leptonic top tagger, but with a limited yield [15] . One strategy to improve the tagging/rejection 
efficiencies is to include the tracker- level mini-isolation observable pr,i/ ^2 Pt, charged > 0.9 summed over all 
charged tracks within a cone around the lepton or muon |16j . 

Our strategy is similar to those in the sense that we attempt to reconstruct (parts of) the top 4-momentum 
without making use of the measured missing transverse energy. This allows us to contrast the top 4-momentum 
with the measured missing energy vector as a background rejection cut. However, we follow a different path to 
extract information on the top 4-momentum. To describe the 3-momenta of the top decay products we define 
the orthonormal coordinate system 

Pbi 
\pm\ 



P° = t~t—\ leading pte direction in b — I decay plane 



HI Pi-iPfP )? 
Wt - (fie ■P D )P B \ 



p" = — — — , m „ ni subleading direction in b — I decay plane 

p ± = p D x p" subleading direction to b — i decay plane. (3) 

Within this coordinate system the neutrino momentum is parameterized as 

Pv = xbP D + x || p 11 + x±p^. (4) 
In the appendix we show the relative size of these three components. 



To reconstruct the neutrino momentum, assuming we know the neutrino mass, we need three components. Be- 
cause we want to compare the neutrino momentum from the reconstructed top to the measured two-dimensional 
missing transverse momentum fir these are three actual unknowns. The W and top mass constraints give us 
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two, so we need to make at least one assumption to describe the collinear decay of a boosted top. As shown in 
the appendix, we can choose two well motivated approximations 

xj_ = decay plane approximation 

xii = orthogonal approximation, (5) 

and correspondingly estimate the neutrino momentum as 

(6) 



Based on this estimated neutrino momentum we can construct measurements which distinguish for example 
the stop signal from the top pair background where the main difference is the existence of additional sources 
of missing momentum in the signal. 



Before constructing such observables we need to sketch a few features of the two approximations for the 
neutrino and top momentum estimates. For example, the mass constraints are quadratic constraints, so we 
cannot expect a unique solution for the neutrino and top momenta. Instead, there will be a set of four discrete 
solutions from which we need to extract the physical solution. The decay plane approximation, which as we 
will discuss in the appendix is motivated by the Jacobian peak in the top momentum we consistently accept 
the solution which returns the smallest top momentum. For the orthogonal approximation, which rejects the 
ti background most efficiently, we take the solution which gives the smaller subtracted missing momentum 

Reconstructing the direction of the neutrino 3-momentum or the top 3-momentum is of course not a full 
reconstruction as we know it from the hadronic top decays using the HEPTopTagger. However, even with 
its reduced numbers of observables it is sufficient for example to reconstruct the azimuthal angle between the 
estimated neutrino or top direction and the measured missing transverse energy 

A</> ee A0(^ T ,p cst ) = 4>(Pr) ~ <Mp cst ) ■ (7) 

The estimated momentum entering this azimuthal difference is not uniquely fixed. We can use the estimated 
neutrino momentum, the estimated leptonically decaying top momentum, or the momentum of the lepton- 
bottom system. 

From the discussion in the appendix estimating the top direction appears to be more promising than es- 
timating the neutrino direction, simply because the neutrino momentum with its large relative uncertainty 
constitutes only a small part of the top momentum. In addition, the studies in the appendix suggest to use 
the orthogonal approximation for this azimuthal angle study. As a starting point we look at the correlations 
between the measured missing transverse energy vector and the reconstructed leptonic top directions. It turns 
out that this observable has the best yield for the stop pair search. In the top row of Figure [T] we show the A<fi 
distributions for the stop signal and the top backgrounds for the orthogonal approximation defined in Eq.([5]). 
It shows a clear difference between signal and tt events for all three choices. The signal shows hardly a peak 
because the reconstructed momentum does not correspond to a physical momentum and is hence uncorrelated. 
The width of the background peak which we need to cut out is smallest when we include the estimated top 
direction, as compared to the neutrino or bottom-lepton directions. This includes a hadron level simulation, but 
the picture might of course change once we include a detector simulation. Qualitatively, all three reconstructed 
directions work at a similar level. 



Improving over this one-dimensional A(f> distribution we can include the correlation with the measured missing 
transverse momentum. It exploits the following feature: for semileptonic top pair events, if the originating top 
is not boosted enough, there exists an upper bound for the neutrino momentum. To observe large ^r, the 
originating top momentum has to be aligned with the ]/>t direction. Again, in Figure [TJ we show the results for 
three reconstructed momentum directions, which work qualitatively the same. 

These observations we can turn into an analysis strategy to extract semileptonic stop pairs after applying the 
five generic steps shown in Table [Ij First of all, it is not guaranteed that given a measured bottom and lepton 
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Figure 1: Top: one-dimensional A(f> distributions with an additional cut p'r > 200 GeV for the stop signal (red) and 
the top pair background (black). Below: two-dimensional p'r vs. A<f) distributions for the it background (2nd row) and 
for the tt* signal (3rd row). From left to right we show A(j>(pr,pi~i), A<l>(pT,Pv) and A.(f>(pT , Pbi) ■ All plots are at the 
hadron level and based on the orthogonal approximation. We apply the universal selection cuts 1-5 and the border of 
the cut 8 is shown in the two-demensional plots. 

momentum configuration we find a valid reconstructed neutrino momentum. As a consistency requirement 
we should check if the reconstructed neutrino momentum combined with the observed bottom and lepton 
momentum really corresponds to the kind of boosted top quarks which we assume to exist in the first place: 

6. possible solution for respective approximation. 

7. pff t > 200 GeV. 

The orthogonal and decay plane approximations perform very differently under the first of these two require- 
ments. By construction and as explained in the appendix, the signal efficiency of the orthogonal approximation 
ranges around 50% while the decay plane approximation almost always delivers a valid solution for the neu- 
trino momentum. However, from Table [TT] we see that this larger signal efficiency is at least compensated by 
the also larger background efficiency, so that in the final result the orthogonal approximation delivers a better 
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Table II: Signal (for different stop masses) and backgrounds after the leptonic top tag. We show both approximations 
to the top momentum defined in Eq.|5|. The j>T vs. A(j> cut is defined in the text. All numbers given in fb. 



signal-to-background ratio as well as a better significance for an integrated luminosity of 20 fb -1 . The situation 
might well be different for smaller luminosities where the number of signal events after all cuts becomes the 
main hurdle. 

As the last step of our analysis, following Figure [T] we require a two-dimensional condition for the azimuthal 
separation defined in Eq.|7|) correlated with the measured missing transverse momentum: 

8. one tagged leptonic top fulfilling all three conditions 

|A *' > I* iOCilev" * T > 200 GeV ^ > T 

All signal and background cross sections after these signal specific cuts we show in Table [TT] As long as we are 
only interested in the direction of the leptonic top momentum but not its size the orthogonal approximation 
works better, which is also shown in Figure [4] in the appendix. Table |H1*| shows that over a wide range of stop 
masses we find of the order of 20 signal events for an integrated luminosity of 20 fb -1 at a collider energy of 
14 TeV. The signal-to-background ratio from the orthogonal approximation ranges around S/B ^1 — 2 and 
the Gaussian significance exceeds S/VB ~ 5 for stop masses from 340 to 540 GeV. For larger stop masses the 
kinematics for the signal as compared to the softer continuum top pair production becomes more distinctive, 



so our analysis should be adapted to this parameter regime. Note that the numbers in Table III are computed 



at the hadron level. Applying all cuts 1.-8. gets rid of essentially all semileptonic top pair events at parton 



level, so the events shown as tt background in Table III are mostly mis-measured leptonic top pairs 



Also in Table |III| we see how the decay plane approximation performs slightly worse for the signal-to- 
background ratio, but gives a similar significance based on a slightly larger number of signal events after 
all cuts. While this significance for semileptonic stop decays might be slightly smaller than for the purely 
hadronic mode [5] the signal-to-background ratio is better and the efficiency of the critical hadronic top tagger 
is reduced to a linear dependence. Combining the two analysis strategies means that based on the order of 
10 fb -1 of LHC data at 14 TeV (or an adjusted number for a lower collider energy) we can discover a top squark 
mainly decaying to top quarks and a dark matter agent in both of the two dominant top decay channels. This 
not only would give us faith in such an observation, it is also the starting point for detailed studies of such a 
new physics signal. 
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Table III: Results based on the reconstructed top momentum direction for the two approximations to the neutrino 
momentum, Eq.j5j. All eight cuts shown in Tables |T| and |ll| are included. The LHC energy is 14 TeV. 
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Figure 2: Mti distributions for the stop signal (red) and the ti background (black). From the left we show the orthogonal 
and the decay plane approximations without the A<f> vs fa cut but with fa > 200 GeV and then both of them with the 
A(f> vs fa cut. The expected endpoint is Mt2 < — m ^.o/ m {) = 418 GeV. 



III. RECONSTRUCTING THE TOP 4-MOMENTUM 



The strategy described in Section [IT] we can also use to fully reconstruct the 4-momentum of the leptonically 
decaying top quark. An additional question becomes how much this additional observable helps in extracting 



the stop signature, as compared to the results shown in Table III Up to now the downside to our stop pair 
analysis is that we do not have a handle on the stop mass. 

If as an additional step we want to gain more information on the heavy particle decaying into a boosted top 
the Mt,2 endpoint is the prime observable and requires full knowledge of the top 4-momentum [55]. As for the 
purely hadronic stop decay the crucial observation is that this stransverse mass endpoint to work we need a 
one-step stop decay into the invisible dark matter agent. A pair of top squarks decaying into tagged tops and 
two dark matter particles is precisely of this topology. 

For the endpoint measurement, we best rely on a combination of the orthogonal and decay plane approxima- 
tions. First, we use the orthogonal approximation for the basic cuts 1.-5. and the more specific top rejection 
cuts 6.-8. Then, we use either the orthogonal or the the decay plane approximation to compute the 771^2 
distribution, i.e. we define a mj 2 variable in terms of the estimated top momenta and the estimated missing 
momentum from other sources than the neutrino as 



T2 



= M T2 ( P : 



HEPTop 
t,h ■ 



Once we measure the leptonic and hadronic top there should appear the endpoint 

M T 2 < mj [ 1 



xi 



(8) 



(9) 



in the stop signal events. In Figure [2j we show the Mj- 2 signal and background resolutions in the orthogonal 
approximation and in the decay plane approximation and clearly see the signal endpoint around 418 GeV. 
From a purely signal point of view and without the two-dimensions Acf) vs. fa cut (8.) the endpoint of the 
stransverse mass defined in the decay plane approximation is more precise. The reason is that in the decay 
plane approximation the estimated top momentum is by construction almost always smaller than the true 
value while in the orthogonal approximation, the estimated top momentum frequently exceeds the true value, 
which spoils the endpoint. Note that the signal-to-background ratio S/B is still fixed by the numbers for the 
orthogonal approximation which we use for the event selection. This signal feature survives after we include 
the final A</> vs. fa top rejection cut. 

The top pair background at hadron level but without the two-dimensional A<j) vs. fa cut is slightly harder for 
the decay plane approximation than for the orthogonal approximation. Again, this feature generally survives 
the application of the A<p vs. fa cut, but with a big caveat arising from the experimental resolution which 
should dominate the behavior of this large- Mx% tail. From the right two panels of Figure [2j we expect that 
the decay plane approximation should work slightly better after applying all base cuts and top rejection cuts, 
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Table IV: Signal (for different stop masses) and backgrounds for the full stop pair analysis. All numbers given in fb. 

but the final word on this comparison needs to include realistic detector effects. At this stage we rely on the 
orthogonal approximation to determine the direction of the top, but on the decay plane approximation for the 
magnitude of the top momentum. 

As the final step of our stop pair analysis we need to check if the Mt2 endpoint is useful for the actual signal 
extraction, i.e. if it improves either the significance or the signal-to-background ratio. The results of such a cut 
we show in Table [TV] We see that the azimuthal separation of the top momentum from the missing energy vector 
alone is highly efficient and gets hardly improved by the additional Mt2 cut. The only reason to consider an 
Mt2 cut in the extraction of the stop signal would be to replace for example the azimuthal direction cut, which 
means replacing the leptonic top reconstruction of the azimuthal angle by the complete top reconstruction by 
our tagging algorithm. Again, this might be an option after taking into account all detector effects, but at the 
level of our analysis the reconstruction of the azimuthal top direction is more efficient. 

IV. OUTLOOK 

In this study we have corrected an earlier statement that at the LHC it is not possible to observe semileptonic 
stop pairs over a wide range of stop masses. In the spirit of the hadronic stop pair analysis [6] we develop a 
search strategy for stop pairs decaying into a boosted semileptonic top pair plus missing energy. We assume 
14 TeV center of mass energy and of the order of 20 fb -1 of integrated luminosity. Different scenarios in collider 
energy and luminosity can easily be explored. 

First of all, the hadronic top decay we fully reconstruct using the standard HEPTopTagger. Earlier studies 
of boosted leptonic top quarks [15|ll6j have shown that when it comes to reconstructing the full top 4-momentum 
the leptonic decay is not too promising. Our approach therefore first focuses on a more limited measurement, 
namely the azimuthal direction of the leptonic top quark. This directional vector we can reconstruct based 
on a decomposition of the neutrino momentum into a leading boost direction and two orthogonal sub-leading 
directions. We show a detailed comparison between two different approximations and come to the conclusion 
that it is promising to reduce the number of kinematic unknowns in the top decay by fixing the neutrino 
momentum orthogonally to the lepton-bottom decay plane. 

Based on this strategy we extract the stop pair signal from the dominant irreducible top pair background 
with 5cr significance for 20 fb -1 of LHC data with a signal-to-background ratio above unity for stop masses 
from 340 to 540 GeV. Our analysis is not optimized for any stop mass, so this result should be very generic 
and should be improved in particular for large stop masses. The discovery potential for semileptonic stop is of 
the same order as for purely hadronic stop decays, both of them based on boosted top quark reconstruction. 

A full reconstruction of the top 4-momentum is considerably harder. On the other hand, to measure the 
stop mass for example using the classical rriT2 endpoint we develop such a strategy. While this additional 
measurement has little impact on the significance or the signal-to-background ratio with which we can extract 
the stop signal it should allow us to determine the stop mass from an niT2 endpoint based on a combination of 
the orthogonal and decay plane approximation of the leptonic top decay. 

The details of the stop mass measurement and its uncertainties need to be determined by an appropriate 
detector simulation. The endpoint value will give us a stop mass measurement in combination with all other 
supersymmetric mass measurements at the LHC. As mentioned in the introduction, the stop mass is particularly 
useful when we combine it with the Higgs mass because the two probe a similar set of model parameters [23] . 

The software used in this study will be published as part of the HEPTopTagger analysis tool [6]. 
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p = 200 - 300 GeV 




p = 400 - 500 GeV 




p = 600 - 700 GeV 



p = 800 - 900 GeV 




Figure 3: (x\\,x±) distributions as defined in Eq.(A3l for top pair production. The four panels correspond to slices in 
the top momentum 200 - 300 GeV, 400 - 500 GeVTBOO - 700 GeV, and 800 - 900 GeV. All events fulfill j> T > 150 GeV. 



Appendix A: Tagging leptonic tops 

In analogy to Ref. [6] we illustrate in detail the motivation of the leptonic top tagger, its implementation, 
and its performance. All results shown in this appendix are based on a hadron level analysis, but without any 
detector simulation, just as the analysis presented in the main body of this paper. In this appendix we focus 
on the reconstruction of semileptonic top pairs, i.e. we only apply the leptonic top tagger to a sample of top 
pair events. 

As mentioned in the main text the W and top mass conditions can be written as 

(Pi +Pb+Pv) 2 =m 2 t 

[pi + Vuf =m 2 w . (Al) 

All three masses of the decay particles are known, which means that once we measure the lepton and b 
4-momenta the allowed neutrino 3-momentum lies on the line which we will show can be described as the 
intersection of an elliptic surface and a paraboloid. 

Following Eq.([3]), repeated here for convenience with pu = pb + pe, 

p B = j-j-j leading pu direction in b — I decay plane 

\PM\ 

p" = 7^ ^r»N^r»i subleading direction in b — i decay plane 

\P£~ (Pe-P D )P D \ 

p ± = p D x pH subleading direction to b — I decay plane. (A2) 

we can write for example the neutrino momentum as 

Pv = xup D + xyp" + x^p 1 - = Puixnp + x\\p" + x^p- 1 ) . (A3) 
In Figure [3] we show the two subleading x distributions for different in the originating top momentum. 

In this parameterization, with = (pe+Pb) 2 , the top mass constraint becomes an ellipse in three dimensions 
(Ptt + m li)x 2 + (pit + fnl t )x\ + m^xn - (m 2 t - ml e )\p be \x D = (m 2 t - m^) 2 /4 

Rh + Rh + Rl, ' (A4) 

with the radii and focal point coordinate 



_ m 2 t - m 2 M _ V\Pbe\ 2 +m 2 M _ \pm\ 

ftt,l — ~ i "-t.2 — *H.l i Xb — -Kt,l ■ (&b) 

2m M m M m u 
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For the W mass constraint it is useful to define the rotated basis vectors 

yu\ _ ( cos6 Mj e sm9 w A (x B 
y\\) \-s\i\Q U j cos 6u,i J \x\\ 

using cos 8i,e. £ — Pe • P° ■ In these variables the W mass constraint reads 



(A6) 



m 4 

\pA 2 V\\ + \Pe\ 2x ± ~ ™w\Pi\VD - U 

.,2 



4 



V\\ , *1 2J/D 



OT TP ^ "nT ft = ^ (A7) 

-"w K w K w 



with the radius Rw — m 2 y/'Api\- 



We can now study the two approximations for the neutrino momentum defined in Eq.([5|. In the decay plane 
approximation xj_ — the two mass constraints become an ellipse and a parabola in the (x\\,xjy) space 

^ + ^_M! = 1 and #-f^ = l. (A8) 

The neutrino momentum we obtain as the intersection of two lines, where in general there exist up to four 
solutions. We use the solution which gives the smallest top momentum. In the orthogonal approximation 
a; || = the mass conditions arc the two ellipses 

gh+ (XD p2 ^ )2 =l and + {X °-f )2 =l, (A9) 

with the rotated radii and focal coordinate 

, _ Rw p// _ Rw , _ cos 8 M/ R 3 

smUbe,e sin ou,e sin Ote,e 

We find these solutions for the orthogonal approximation numerically, being aware of the fact that in some cases 
there does not actually exist any solution. This finite efficiency we observe in our analysis results in Table \U\ 
In case we find more than one solution we keep the one which gives the smaller subtracted missing momentum 

h-p c T % 

In Figure [4] we show several measures describing the difference between the estimated leptonic top and 
neutrino momenta and the respective true momenta at parton level: 



A os V,| 


= Ku\ - 


\p17 




= <fC - ' 




A CBt e t . u 


= c - ( 


itruc 



(All) 

At hadron level the orthogonal approximation shown in red performs consistently better than the decay plane 
approximation to reconstruct the direction of the neutrino and the top momenta. We also see that the azimuthal 
direction of the top quark can be more precisely reconstructed than the neutrino's azimuthal direction. This 
occurs because compared to the bottom and lepton reconstruction the neutrino is the weak spot of the top 
reconstruction. This is why in our analysis described in Section|n]we use the top azimuthal direction estimated 
in the orthogonal approximation. Only considering tops with px.t > 300 GeV only slightly improves this 
picture. 

If we do not apply the condition j>T > 150 GeV effectively on the neutrino momentum the corresponding 
analysis to Figure [4] paints a different picture: the decay plane approximation achieves the better resolution at 
parton level, and this features is largely the same at hadron level. This is the reason why our Mt.2 endpoint 
analysis is based on a mass reconstruction in the decay approximation (after the event selection in the orthogonal 
approximation) . 
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Figure 4: Absolute and relative A est |p|, A est |p|/|p truc |, A cst <j) and A cst 8 distributions in the decay plane approximation 
(black) and the orthogonal approximation (red). The upper two rows show the neutrino momentum estimate, the lower 
two rows the top momentum estimate. The first of the two rows is at parton level, the second row at hadron level. Only 
semileptonic top pair events are included. At hadron level we only require j>T > 150 GeV while at parton level all events 
also fulfill p T ,i > 20 GeV and p T ,b > 25 GeV. 



Appendix B: Jacobian in the decay plane approximation 



The relative position of the top decay products is completely determined by the vector of variables X = 
(7t, cos 6b, cos 8,(f>). Here, -f t is the boost factor of the top in the laboratory frame p\ ah = /3 t jtmt, Ob is the 
t — > bW decay angle with respect to the top direction, and (9, 4>) are the W decay angles in the W rest frame. 
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In terms of these variables, the three momenta of the decay products in the top rest frame are given by 



Pb 




Pi 



m w 



H — 7/3 cos 9\ 
sin 9 cos <j 
sin 9 sin q 
\7 cosQ — j(3j 



Pv 



m w 



/ 7 + 7/3 cos 9 

— sin 9 cos 4> 

— sin sin <fr 
\— j cos 9 — 7/3/ 



(Bl) 



P and 7 are defined by the W and 6 momenta in the top rest frame 



fiimw = \Pw\ = \Pb\ 



m t 



'1- 



1 - 



m t 



(B2) 



and -Eb = \fp^+ m\. For an unpolarized top the top decay phase space distribution in terms of the variable 
vector X can be written as HU 



P(X) = P t ( 7t ) iV(cos< 



(B3) 



where P t ~ exp(— const x j t )/j t is the distribution for the top momentum and Pw is the distribution for the 
W decay angle 9. 

The problem of P(X) is that it does not correspond to our measurements listed in Eq.([2]). Linking the 
observables X to our energy measurements in the laboratory frame returns a Jacobian which dominates the 
distribution of the top decay products in our analysis. After boosting the lepton-bottom system by j t in the 
direction of 9b the lepton and b energies and their invariant mass from Eq.Q become 

plab 



Ef a = lt E b + lt p t \p b \cos9 b , 

7i7(l — cos 9) ~ jtf3 t sin 9 b sin0sin</> + ^ t fitl cos9 b (cos9 — /?) 



^lab = m W 



mf — m w + m 



m t \p b \ cos 9 , 



(B4) 



These observables we can express in terms of dimensionless vector Y — (7^, E^ )1 Eg, Th^g) including 



77'iab 

fp _ b 

\Pb\ 



E t = 2- 



mw ' 



'be 



'hi 



m t \pb\ 



(B5) 



Using the relation P(X)dX — Py(Y)dY, the distribution Py(Y) after a variable transformation can be calcu- 
lated from the identity 



P Y (Y) =det 
with det 



dY 



dX 
dY 



<-w 



dX 



P < MP «\ »*|ft| 

7t 2 01 sin# b sin 9 | cos0| = Tt( S t y %f + 2B lt + C 



in terms of 9 given by the last line of Eq. ( B4 ) and 

A = - sin 2 



El 



\Pb? 



-1 - 7 2 sin 4 0(l + /3)^ 



B = sin 2 



E b E b 

IpW 



- 7 sin 2 9(1 + 13) [E e - 7 (cos 9 - f3)E b ] 



C = - sin 2 9(1 + El) - [Ei - 7(cos6> - f3)E b } 2 , 



(B6) 



(B7) 



Once we measure the Y components (E b , Eg, rhbi) we can estimate the top momentum according to this proba- 
bility Py(Y). The poles in the Jacobian det \OY/dX\ = give us a pronounced Jacobian peak and constrain the 
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allowed range of the top momentum. Following Eq.(B6) the peak position implies cos0 = 0, which corresponds 
to the kinematic configuration where pb, pi and p v are coplanar. 

Moreover, the pole condition is a quadratic equation for y t which has two solutions. Because Pt{lt) is steeply 
falling for increasing j t the most likely value for "f t is the smaller one. To estimate the top quark momentum we 
therefore choose the solution which yields the smaller momentum in the decay plane approximation (xj_ = 0). 
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